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Abstract : Numerical integration is a frequently-needed tool 
in modern Science and Engineering. Engineers and Scientists 
typically visualize integration as the process of determining the 
area under a curve. Besides this, because of its many more 
applications, it is often viewed as a discipline in and of itself. 
In this paper we develop a mathematical simulator for solving 
numerical integration problems. This simulator is incorporated 

with a combination of Trapezoidal rule, Simpson's 1/3 rule, 

Simpson's 3/8 rule and rule. 
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1. Introduction 

Integration generally means combining parts so that they form a 
whole. The foundation for the discovery of the integral was laid 
by Cavelieri, an Italian mathematician, in around 1635. Besides, 
numerical algorithms are almost as old as human civilization. In 
numerical analysis, numerical integration constitutes a broad 
family of algorithms for calculating the numerical value of a 
definite integral. Archimedes of Syracuse (287-212 BC) created 
much new mathematics, including the "method of exhaustion" 
for calculating lengths, areas, and volumes of geometric figures 
[ I When used as a method to find approximations, it is in much 
the spirit of modern numerical integration; and it was an 
important precursor to the development of the calculus by Isaac 
Newton and Gottfried Leibnitz. Following Newton, many of the 
giants of mathematics of the 18 th and 19 th centuries made major 
contributions to the numerical solution of mathematical 
problems . 

Given a set of data of points (x ,y ), (x { , y { ) , (x i ,,y il ) 
of a function y = f( x ) where f(x) is not known explicitly, 
it is required to compute the value of the definite integral 



\ydx 



(1) 



We derive a general formula for numerical integration using 
Newton's forward difference formula. 

Let the interval \a,b\ be divided into n equal subintervals such 
that a - X < X r < ... < X n = b. 



Clearly, x n = x Q +nh ■ Hence the integral becomes 

I = j y dx ■ No w approximating y by Newton's forward 

difference formula, we obtain 131 

, fr A P(P _ 1) a2 P(P- !)(/>- 3) .3 i , 

I=\[y + pAy + lu J A 2 y + FK1 -ff -A y + ...] dx 

J 1.2 l.Z.j 

Since X = X + ph and so dx = dp hence the above integral 
becomes 





This implies 



1.2 



1.2.3 



\y dx = nh[y +(n/2)Ay + — A 2 y + ... ]dp 

J 12 



(2) 



From this general formula, we can obtain different integration 
formula by putting n = 1, 2, 3, .. . etc. 

Programming language C is very flexible and powerful. It 
originally designed in the early 1970s . It allows us to 
maximum control with minimum command. It is recognized 
worldwide and used in a multitude of applications especially in 
numerical analysis. Along with other numerous benefits, we 
have used programming language C in this paper. 

The outline of this paper is as follows: Section 2 contains the 
brief description of the existing methods with methodology and 
simple examples. In Section 3, we develop a numerical 
simulator, using the programming language C, which gives us 
the solution of a problem simultaneously regarding four popular 
existing numerical integration methods namely Trapezoidal rule, 

Simpson's 1/3 rule, Simpson's 3/8 rule and Weddle's rule. 
Moreover, the simulator identifies the method that gives the best 
solution comparing with possible exact solution of the problem. 
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We devote Section 4 to an output of the program for a specific 
problem. Conclusions are given at the end at Section 5. 

2. Existing Methods 

We give a brief description of the existing methods of numerical 
integration like Trapezoidal rule, Simpson's 1/3 rule, 
Simpson's 3/8 rule and Weddle's rule in this section with their 
methodology and simple examples. 

2.1 Trapezoidal Rule 

Putting n—1 in (2) all differences higher than the first will 
become zero and therefore we have 



■M 

Similarly, | y dx = (h/3)[y 4 + 4y 3 + y 2 ] and finally 
\y dx = (h/3)[y n +4y n _ 1 + y n _ 2 ] 

x n-2 

Combining all these expressions, we obtain 

X n 

\y dx = (h/3)[y Q +4(y^ + y^ + ■■■+ y„ _^) + 



1 ' J 3 



y n-Y 



■*] 

Jy dx = h[y + (l/Z)Ay \ = (h/2)[y + y 1 \ 



+ 2(y 2 + y 4+ ... + y n _ 2 ) + yJ 

This is known as Simpson's 1/3 rule. 
For the next interval [x,, x 2 ] and others we have the similar 2.2.1 Example 

expression as j" y dx — {hi 2) [y x + y 2 ] an< ^ so on an ^ fi nairv 
for the last interval [■*„_] , X n ] , we have 



Calculate the integral of the function, f(x) — 2x in the 
interval (0,1) using Simpson's l/3rule. 

Let h = 1/ 3 . The tabulated values of X and y are given below 



\ydx = (h/2)[y n _ 1 + y n ]. 

Combining all these expressions, we obtain the rule r21 
j y dx = (h 1 2) [ y Q + 2( y, + y 2 + . . . + y„_, ) + y J 
This is known as Trapezoidal rule. 
2.1.1 Example 



X 





1/6 


1/3 


Vi 


2/3 


5/6 


1 


y = f(x) 





1/3 


2/3 


1 


4/3 


5/3 


2 



Using Simpson's 1/3 Rule 



1 



+ 4 



1 



2 4 



We are estimatin 



g f 1 dx using the Trapezoidal Rule 

J , lnx+3 



with n — 4. 

Here h = (b - d)ln = (3 + l)/4 = 1, and 

f- -^ = (/i/2)[y +2(y 1 + y 2 +...+ y„_ 1 ) + y„] 

J j In x + 3 

= [0.9102 + 2(0.721 + 0.621 + 0.558) + 0.514]/ 2 = 2.613 
2.2 Simpson's 1/3 Rule 

Putting Tl — 2 in (2) all differences higher than the second term 
will become zero and therefore we have 

x 2 

jydx = 2h[y +Ay + A 2 y,/6] = (h/3)[y 2 + 4y + y ] 



2.3 Simpson's 3/8 Rule 

Putting n — 3 in (2) all differences higher than the third term 
will become zero and therefore we have 



3 3 9 9 3 3 

J y dx = 3/z[y„ + — Ay„ H A y„ H A y n ] 

L -^0 2 12 24 



= -^[^0+3^+3^ + ^] 

Similarly, f y dx — —h [y 3 +3y 4 + 3y 5 + y 6 ] and so on. 
J 8 

Summing up all these, we obtain f3 ' 

J j <ij<; = (3/z/8)[j + 3^ +3y 2 +2y 3 +3y 4 +3y 5 + 

+ 2 y b + ... + 2 j„_ 3 + 3 y n _ 2 + 3y„_! + y„ ] 
This is known as Simpson's 3/8 rule. 
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2.3.1 Example 

We solve the integral of the function, f(x) = x 2 +1 , in the 
interval [1,4] using Simpson's 3/8 rule. 

Let h — 1- The tabulated values of X and y are given below: 
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X 


i 


2 


3 


4 


y 


2 


5 


10 


17 



From Simpson's 3/8 Rule, 

jy dx = (3h/S)[y + 3y l + 3y 2 + 2y 3 + 3y 4 + 3y 5 + 

+ 2y 6 + ... + 2y n _ 3 + 3y n _ 2 + 3y n _ l + yj 
2.4 Weddle's Rule 

Putting n — 6 in (2) all differences higher than the first will 
become zero and we obtain 

\y ^=7^[3 ; +5 > 7 l + >'2 +6 >'3 + >'4 +5 >'5 + >'6] 
J 1(J 

And thus the general form for Weddle's rule is 

x n 

| y dx = (3h/10)[y +2(y 6 + y n +...+ y n _ 6 ) + 
x 

+ 5(y l + y 1 + ...+ y) + 6(y 3 + y 9 + ... + y n _$) + y n ) 
2.4.1 Example 

6 

We evaluate ^2x dx U P to three decimal places, by Weddle's 

o 

rule. 

Let h — 1.0. The tabulated values of X and y are given below 



X 





i 


2 


3 


4 


5 


6 


y 





2 


4 


6 


8 


10 


12 



From Weddle's Rule, 

x„ 

jy dx = (3h/l0)[y +5(y 1 + y 5 ) + 6y 3 + 

+ (y 2 + ^4) + ^] = 36 
3. Algorithm 

Now we are introducing an algorithm to evaluate an integral 
numerically by four well known methods namely Trapezoidal 
rule, Simpson's 1/3 rule, Simpson's 3/8 rule and Weddle's rule 
simultaneously and it will ensure which method gives the best 
accuracy compare to the exact solution. 
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INPUT: function f , limits X Q and X n , number of division n , 
direct result r . 



Step-1: Compute h ■ 



X n X () 



Step -2: Set i - 

Step-3: While i<n, repeat Step- 4 
Step-4: Set y [i] = f(x + ih) 
Step-5: Set i - 1 

Step-6: While i < ft, repeat Step-7 
Step-7: Set sum _1 = sum _1 + 2y [i] 
If i % 2 = 

Set sum _ 2 = sum _2 + 2y [i] 
Else 

Set sum _ 2 = sum _2 + 4y [/] 
If z%3 = 

Set sum _ 3 = sum _3 + 2y [i] 
Else 

Set sum _ 3 = sum _ 3 + 3y [i] 
t_Sum_l = (h/2)*(y[0] + y[n] + sum_l) 
t_Sum_2 = (h/3)*(y[0] + y[n] + sum_2) 
t_Sum_3 = (3h/S)*(y[0] + y[n] + sum_3) 
If i%6 = 

Set sum _ 4 = sum _4 + 2y [i] 
Else if i % 3 = 
Set sum _ 4 = sum _ 4 + 6y [i] 
Else 

Set sum _ 4 = sum _ 4 + Ay [1] 
Else If (i%6 = llli'%6 = 5) 
5«m _ 4 = 5«m _ 4 + 5y [i] 
Else 

5«m _ 4 = sum _4 + y [i] 

t_sum_4 = (3*h/8)*(y[0] + y[n] + sum_4) 
Step-8: Set a — \t _sum _\ — r\ , b — \t _sum _2 — r\, 
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c = \t _ sum _ 3 — r , d = f _ .swm _ 4 - r 

If ( a < b and a < c ) 
if (a<d) 
? _ Sum _ 1 
Else 

t _ sum _ 4 

If ( b < a and b < C ) 

If (&<d) 

£_Swm_2 
Else 

? _ sum _ 4 

If ( C < a and C <b) 

U(c<d) 
t _Sum_3 
Else 

t _ sum _ 4 

OUTPUT: t_Sum_\, t_Sum_2, t_Sum_3, 
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t _ sum 
accurate. 



4 with message which sum is most 



STOP. 



4 Findings 



We are now giving a trial of the simulator for the numerical 
integration problem / = f dx ■ 

il + x 



4.1 Input: 

This is a program to evaluate numerical integration and check 
which method gives the best accuracy. 
Enter equation:- 
1 



1 + x 

Type the initial value & final value: 



Please enter 6 or multiple of 6: 
12 

Type direct result of integration: 

.693140 

4.2 Output: 

The Trapezoidal rule gives: 0.693581 
The Simpson l/3rule gives: 0.693149 
The Simpson 3/8 rule gives: 0.693150 
The Weddle's rule gives: 0.693147 
Weddle's rule gives the best among these. 

5 Conclusion 

In this paper, we develop a simulator incorporated with 
traditional Trapezoidal rule, Simpson's 1/3 rule, Simpson's 
3/8 rule and Weddle's rule for solving numerical integration. 
We observed that the result obtained according to our procedure 
is completely identical with the hand calculation. We therefore, 
hope that this simulator can evaluate integrals numerically to 
save time and labour. Moreover, we have seen that Weddle's 
rule gives the best solution among the four methods. 
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Type the number of division. 
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